New Algorithm for Parallel Laplacian Growth by Iterated Conformal Maps 
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We report a new algorithm to generate Laplacian Growth Patterns using iterated conformal 
maps. The difficulty of growing a complete layer with local width proportional to the gradient of 
the Laplacian field is overcome. The resulting growth patterns are compared to those obtained by 
the best algorithms of direct numerical solutions. The fractal dimension of the patterns is discussed. 



Laplacian Growth Patterns are obtained when the 
boundary T of a 2-dimcnsional domain is grown at a rate 
proportional to the gradient of a Laplacian field P pj. 
The classic examples of such patterns appear in viscous 
fingering in constrained geometries (like Hele-Shaw cells 
or porous media). Here a less viscous fluid (inside the do- 
main bounded by T) displaces a more viscous fluid which 
is outside the domain. The field P is the pressure, and 
Darcy's law determines the velocity v to be proportional 
to VP. Using the incomprcssibility constraint outside 
the domain V 2 P = 0, and each point of T is advanced at 
a rate proportional to VP |2|, |3j ■ Thus in numerical algo- 
rithms 0| one needs at each time step to add on a whole 
layer to the pattern with a local width proportional to 
|VP|. The boundary conditions are such that in radial 
geometry as r — > oo the flux is VP = const x f/r. On the 
boundary T one usually solves the problem with the con- 
dition P — UK where a is the surface tension and k the 
local curvature of T Q . Without this (or some other) ul- 
traviolet regularization Laplacian Growth reaches a sin- 
gularity (in the form of a cusp) in finite time In this 
Letter we present a new algorithm to grow such patterns 
using iterated conformal maps. The basic method was 
introduced in the context of Diffusion Limited Ag- 
gregation (DLA) where it was successfully employed 
to solve a number of outstanding problems in the the- 
ory of DLA 0, , its fractal dimension jToL [TT| and its 
multifractal properties 0, 0, 0] . 

At the heart of the method stands the elementary map 
<j)\ e which transforms the unit circle to a circle with a 



"bump" of linear size vA around the point w 
employ the elementary map Q 
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This map grows a semi-circular bump with two branch 
points at the angular positions ± a, where 



a = tan 



1 - A 



(2) 



By iterating this fundamental map with randomly chosen 
angles n and the bump sizes A n chosen such as to obtain 



equal size particles on the cluster, 



A„ = 



$(n-l)' ( e i0„)| 



one can easily grow a DLA pattern. On the other hand, 
the direct application to the closely related problem of 
viscous fingering remained unaccomplished due to tech- 
nical difficulties that have been now surmounted, as we 
report below. 

We are interested in <&( n ^(w) which conformally maps 
the exterior of the unit circle e in the mathematical 
u;-plane onto the complement of the Laplacian pattern 
in the physical z-plane. As in previous work the map 
$( n )(w) is obtained by iteration of fundamental maps 
4>Xj.9j{w)- The superscript n denotes n growth steps. 
The gradient of the Laplacian field VP(z(s)) is 



\VP{z{s))\ = 



1 



z{s) = $(™)(e iy ) . (3) 
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Here s is an arc-length parametrization of the bound- 
ary. Contrary to DLA which is grown serial, i.e. particle 
by particle, in Laplacian growth for the problem of vis- 
cous fingering we need to grow in parallel, i.e. layer by 
layer. Nevertheless the map & n \w) is still constructed 
recursively. Suppose that we completed the last layer at 
growth step m, i.e. $( m )(«j) is known, and we want to 
find the map $( m +p) [w) which maps the exterior of the 
unit circle to the exterior of the pattern after the addi- 
tion of one more layer whose local width is proportional 
to |VP(z(s))|. The number p of growth events should 
be arranged precisely such as to add the aforementioned 
layer. To grow a full layer of non-overlapping bumps us- 
ing the elementary map Q is difficult because the series 
of sizes {A m +/c}^ =1 and of positions {0 m +k}^ = i depend 
on each other. 

Consider then the (rn + k)th step of growth that 
we implement at the angle m +k- Due to the re- 
paramterization during the p growth steps k = 1,2, . . .p 
we need first to find an angle 9 m +k according to the fol- 
lowing rule. For a given position, the size has to be such 
that its image under the map $( m+fe " 1 ) is proportional 



to the local field 



i.e. 
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here 9 m+ k is defined through 



$( m ) 



^(m+fc-l) ^gie m+fc ^ 



(4) 



On the other hand, the position has to be such that the 
new bump on the cluster precisely touches the previously 
grown one, i.e. the image of one of its branch points has 
to be equal to the image of the corresponding branch 
point of the previous bump. When growing the layer in 
a mathematically positive direction this means that 

q>[rn+k)(j,{9 m+h -oi m+h \ _ ^(m+H)^(8 mtl -i+a m+ n)| 

(5) 

where according to Eq.JSJ a n = tan 1 (2\A^/(1 — A„)j 
Note that a n only depends on the size A„ and not on 
the position, because it is a property of the un-rotated 
fundamental map, namely its branch points. By use of 
the inverse fundamental map equation JSJ simplifies to 



t+k . e i(S m + k -a m+k ) 



— e i(e m +fc-i+Q m +fc-i) 



(0) 



where 4>\ m+k ,o is the fundamental map which depends 
only on the size of the bump y/\ m +k- This can be further 
simplified to get an equation for 9 m +k 



oi(6 m +h-X — <Xm+h-l) 

<px m+k ,o pWj 



(7) 



where the right hand side depends only on the previous 
bump (m + k — 1) and the size of the new particle (m + k). 
Thus we can compute 9 m+ k from equation (JJJ once we 
know the size of the bump. 

To accomplish the goal of growing a full layer we now 
need to estimate the size of the bump X m +k and from 
there compute the position Q m +k according to equation 
Q. This will exactly close the gap between successive 
bumps in the layer. In estimating X m +k we use the fact 
that the bumps are minute on the scale of the growth 
pattern. Therefore the field does not vary significantly 
from bump to bump. Thus we set the value of X m +k to 



A 



An 



,%0„ 



(8) 

In other words, we estimate \ m +k for the unknown posi- 
tion 9 m+ k at the position of the last bump 9 m+k _i. Note 
that with equation (JJJ) we are guaranteed to have a full 
cover of the layer independently of the estimate of \ m +k- 
Only the very last bump may not "fit" in, and is therefore 
not grown. 

In order to avoid a bias in the growth pattern we al- 
ternate the direction of filling the layer from clockwise 
to anti-clockwise and vice versa. Furthermore we choose 
the position of the first bump of a given layer randomly 
on the unit circle. Thirdly, to avoid growing very many 
very small particles in the fjords, we introduce a cut-off 




FIG. 1: Viscous fingering pattern obtained through 100 
000 growth events. The intermediate steps of growth are 
also shown. Every drawn intermediate pattern is an actual 
conformal map of the unit circle, with the last one being 

<j>(100 OOOWgtflx 



A cu t for A„. That means that if a bump is about to be 
grown with a smaller A n then we just avoid growing it 
and proceed to a position which is A cu t further down the 
unit circle. The results of our algorithm were checked to 
be independent of A cu t in the range 1CP 6 < A cu t < 10 -10 . 
The results shown in this Letter are for the lowest cut-off 
Xcut= 10" 10 . 

Finally, the conformal map at the end of the growth of 
the layer can be written as 



■°0A m+p m+p M , (9) 



where o stands for a functional composition. In Fig. ^ 
we present 3>( 100 000 \lo) which is a growth pattern with 
100 000 growth events, including the intermediate stages 
of growth. The growth patterns is very similar to those 
obtained by direct numerical simulations of viscous fin- 
gering in a radial geometry. For comparison we show in 
Fig. |21 the pattern obtained by direct numerical solution 
□ 

One advantage of the present approach is that the con- 
formal map 

$(») (u) is given explicitly in terms of an iter- 
ation of analytically known fundamental maps, providing 
us with an analytic control on the grown clusters. The 
map (uj) admits a Laurent expansion 



<i>( n V) 



F 



(>') 



(10) 



The coefficient of the linear term is the Laplace radius. 
Defining the radius of the minimal circle that contains 



3 




FIG. 2: Viscous fingering pattern obtained by direct numeri- 
cal solution, see Q. 



the growth pattern by R n , one has the rigorous result 
that 

Rn > F[ n) > R n /4 . (11) 

It is therefore natural to define the fractal dimension of 
the cluster by the radius-area scaling relation 

F (n) _ S i/D f (12) 
where is the area of the cluster, 

n 

5 (n) = J2 X i ^^"^'(e^')! 2 ■ (13) 

3=1 
(n) 

On the other hand F± is given analytically by 

n 

^"^nvti+^y, (w) 

and therefore can be determined very accurately. In Fig. 
|3]we plot the area in double logarithmic plot against 
Ff 1 ' . The local slope is the apparent dimension at that 
cluster size, it appears that the dimension asymptotes to 
the value indicated by the straight line which is D = 1.7. 

A few comments are in order. First, our algorithm 
does not employ surface tension; rather, we have a typ- 
ical length scale Ao that removes the putative singulari- 
ties. The boundary conditions for the Laplacian field are 
still zero on the cluster. The apparent correspondence of 
our patterns with those grown with surface tension regu- 
larization points out in favor of universality with respect 




J I I I I l_l I I L 



1 10 




FIG. 3: Double logarithmic plot of the area S (n) vs. f\ ( 



to the method of ultra-violet regularization. Second, one 
may worry that growing semi-circular bumps means that 
the elementary step has both length and height propor- 
tional to |VP|, and not just height. Since our elemen- 
tary growth events are so small in spatial extent this does 
not app ear to be a real worry. Lastly, in previous work 
[Pa. Ha. flT| the technical difficulty of growing a complete 
layer was circumvented by constructing a family of mod- 
els which included Laplacian growth only as a limiting 
model. This family of models achieved a partial coverage 
C < 0.65 of every layer of growth, with C = 1 being Par- 
allel Laplacian growth which could be considered only as 
an extrapolation of lower values of C. Assuming mono- 
tonicity of geometric properties as a function of C, the 
extrapolation procedure indicated strongly that viscous 
fingering were not in the same universality class as DLA, 
having an extrapolated dimension D = 2. The present 
algorithm which achieves directly the limit C = 1 ap- 
pears at odds with the extrapolation procedure. The di- 
rect measurement of the dimension is in close correspon- 
dence with DLA whose dimension is 1.713±0.03 At 
present it is not clear whether the contradiction is due 
to a non-monotonicity as a function of C, Whether the 
clusters have not reached their asymptotic properties or 
whether there is another reason that will be illuminated 
by further research. 

At any rate it is hoped that the availability of a di- 
rect procedure to grow viscous fingering patterns by it- 
erated conformal maps will help to achieve a theoretical 
understanding of this problem on the same level of the 
understanding of DLA. 
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